The pair potential of colloidal stars 
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We report on the construction of colloidal stars: 1 polystyrene beads grafted with a dense brush of 1 
long and 10 nm wide semi-flexible filamentous viruses. The pair interaction potentials of colloidal stars 
are measured using an experimental implementation of umbrella sampling, a technique originally developed 
in computer simulations in order to probe rare events. The influence of ionic strength and grafting density 
on the interaction is measured. Good agreements are found between the measured interactions and theoretical 
predictions based upon the osmotic pressure of counterions. 



Polyelectrolyte brushes of flexible polymers have been the 
subject of many theoretical fl','2','3'] and experimental |3] stud- 
ies. Recently focus has shifted to semiflexible brushes |4] 
for which the persistence length P is large compared to the 
monomer separation, but small compared to their contour 
length L, or P << L. In contrast, here we investigate brushes 
with P ^ L. The grafted brushes consist of bacteriophage 
Ml 3 viruses, which are rodlike, semiflexible charged poly- 
mers of length L ~ 880 nm, diameter D — 6.6 nm, and 
persistence length ^ 2fim f3\. The bare, linear charge density 
of M13 is high; ~ 7e"/nm. 

In this letter, we describe "colloidal stars", which are anal- 
ogous to star polymers filSSQl, constructed by grafting ge- 
netically engineered M13 viruses JH] to polystyrene spheres. 
These stiff brushes represent a new class of stars. The M13 are 
rigid enough to form liquid crystals l'9'l, but when grafted to 
a sphere remain flexible enough to be distorted by the direc- 
tor field, as shown in Fig. [T] These colloidal structures are 
characterized by fluorescent microscopy, transmission elec- 
tron microscopy (TEM) and fluorometry. The interaction po- 
tential is probed using laser tweezers. To extract the steeply 
varying pair-potential we develop a new experimental proto- 
col based on the computer simulation method known as um- 
brella sampling [10], but modified to increase the protocol's 
efficiency under experimental constraints. This new method 
allows the measurement of potentials much greater in magni- 
tude than done perviously with line traps flU [12I1 . We find 
that the measured potential of the colloidal stars can be mod- 
eled as arising from the osmotic pressure of the counter-ions, 
which is in several fold excess of the repulsion due to rod ex- 
cluded volume. 

Ml 3 bacteriophage was grown and purified as described 
elsewhere [13]. The M13 capsid protein pIII, present only on 
one end of the virus, was modified in order to display cysteine 
residues. We achieved this through making use of the Ph.D.- 
C7C Phage Display Peptide Library (M13-C7C, New England 
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Biolabs, Beverly, MA) j^. We were able to create colloidal 
stars with core sizes varying from 10 nm to 1 fim. Fig.[T^ and 
b show M13-C7C viruses conjugated with 10-nm colloidal Au 
particles (Ted Pella, Redding, CA). In this article, we focus 
on the colloidal star constructed by attaching the engineered 
phages to a 1 micron diameter polystyrene sphere. This was 
done using the following procedure: First, 230 pd of 8.8 
mg/ml M13-C7C was reduced with 2 ^1 of 0. 18 mg/ml TCEP 
(Tris(2-carboxyethyl)phosphine) for 15 min. This M13-C7C 
solution was mixed with 2 /il of 19 mM maleimide-PE02- 
biotin (Pierce, Rockford, IL) for 1 h in 20 mM phosphate 
buffer at pH = 7.0. The phage solution was dialyzed ex- 
tensively against phosphate buffer to remove excess biotin 
and the pH was readjusted to 8.0. Subsequently, the phages 
were mixed for 1 h with 1 mg/ml Alexa Fluor® 488 car- 
boxylic acid succinimidyl ester (Molecular Probes, Eugene, 
OR), and centrifuged four times at 170, OOOg for I h to remove 
free dye molecules. 0.5 mg/ml of the fluorescently-labelled 
viruses were then incubated with 0.5%(w/v) straptavidin- 
coated polystyrene beads of diameter d = 0.97 ± 0.02 /im 
(Bangs Laboratories, Fishers, IN) for 24 hours at room tem- 
perature. To the suspension 0.05 mg/ml a-casein (Sigma, St. 
Louis, MO) was added and the whole mixture was centrifuged 
twice at 20,000(7 for 10 min. Finally, the pellet was resus- 
pended in 100 fil of phosphate buffer (5 mM, pH 8.0) and 
stored at 4°C. The number of the sphere-bound viruses was 
determined using a fluorescence spectrophotometer (F-2000, 
Hitachi, Tokyo, Japan). By varying the stoichiometric ratio of 
biotinylated viruses to straptavidin-coated beads we created 
star polymers of different grafting densities as revealed by 
both fluorescence and TEM images (Fig.fHc-e)). Fluorescent 
images were taken on a fluorescence microscope (TE2000-U, 
Nikon) equipped with a lOOx oil-immersion objective and a 
cooled CCD camera (CoolSnap HQ, Roper Scientific). The 
TEM samples, stained with 2% uranyl acetate, were imaged 
with a 268 microscope (Morgagni, FEI Company, Hillsboro, 
OR), operating at 80 kV. 

At the grafting density of 135 phages/bead (Fig.[T^), the an- 
chored dye-labelled rods form a spherically symmetric corona 
around the bead with a radially-averaged intensity (RAI) pro- 
file shown in Fig. [TJf)- We model the phage-grafted bead as 
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FIG. 1: (a) and (b): TEM images of 10 nm Au-bound M13 viruses 
of different nanoarchitectures. (c)-(e): TEM (right panel) and fluo- 
rescence (left panel) images of labeled phage grafted to unlabeled 1 
fim PS beads with varying grafting densities, (c) 3 phages/bead. (d) 
38 phages/bead. (e) 135 phages/bead. (f): Radially-averaged fluo- 
rescent intensity profiles of the phage-grafted bead. Symbols: exper- 
iment; Solid curve: theoretical calculation with varying orientational 
order parameters (S) of anchored rods, (g): Fluorescent image of 
colloidal star in a M13 nematic (in contrast to (e) where the solvent 
is isotropic). The "hair" grafted to the bead is "combed" parallel to 
the director by the nematic. (h): Brightfield image of colloidal stars 
associating end-to-end in a M13 nematic. (i): Fluorescent image of 
(h). The combed stars associate in chains aligned parallel to the ne- 
matic director with surfaces separated by a micron. Bare spheres in a 
nematic also assemble into chains, but with surfaces in contact. The 
scale bars are 500 nm. 



hard rods anchored to the sphere with a Gaussian angle dis- 
tribution, which is centered around the surface normal. The 
diffraction-limited fluorescent image of the colloidal star was 
computed by convolving the distribution of the rod's segments 
with the theoretical 3D point spread function (PSF) of the mi- 
croscope |14J. As can be seen from Fig. [Hf) the calculated 
RAI profiles are insensitive to the orientational order parame- 
ter of the anchored rods S ~ ^{3 cos^ 6—1), where 9 is the 
angle between the rod and the surface normal. However, the 
best fits were for intermediate order parameters. 

The free energy as a function of separation between two 
colloidal particles Wint(7') (the potential of mean force) can 



be determined up to an additive offset by the Boltzmann rela- 
tion, P{r) ^ exp[— Wint(r')/fcBT]- Experimentally this is ac- 
complished by measuring the probability P{r) of finding the 
particles at a separation r. However, for states of even mod- 
erate repulsive interaction energies P{r) becomes very small. 
As a result, infrequent visitation of improbable states leads to 
poor statistics and errors in the determination of P{r) which 
limited the magnitude of measured potentials in previous im- 
plementations of line traps, or single bias potentials to about 6 
ksT lfTll[l2ll . In this paper the maximum measured potential 
is 40 feflT, but we estimate that potentials several times this 
value are feasible with the laser power and optical resolution 
of our instrument. 

We achieve these measurements by employing the method 
of umbrella sampling, in which a biasing force is used to 
enhance sampling of rare configurations; results are then re- 
weighted to obtain the physical probability distribution IToll . 
Specifically, we place two colloidal stars (Fig. ^) in sepa- 
rate laser traps and measure the histogram of separation dis- 
tances between the colloids. The measurement is performed 
in a series of windows, each of which uses a different sepa- 
ration distance between the minima of the two laser traps. In 
each window the stars fluctuate about the minimum of a total 
potential resulting from a combination of the dual traps and 
interparticle star potential. Only 6 k^T of each of the total po- 
tentials is sampled and each minimum has a different energy, 
but here we show how the total potentials from overlapping 
windows can be combined to produce a single interparticle 
pair-potential of large range and magnitude. For the protocol 
typically used in simulations, results from different windows 
would be simultaneously re-weighted and stitched together to 
obtain a continuous function for the probability P(r) using the 
weighted histogram analysis method (WHAM) lfl5ll . How- 
ever, the biasing potential is a function of two coordinates be- 
cause the position of each bead is controlled by a separate 
trap. The number of independent measurements required for 
a particular level of statistical accuracy using WHAM rises 
exponentially with the number of dimensions of the biasing 
potential (even if the probability is projected onto a single co- 
ordinate). We overcome this limitation as follows. 

Our goal is to measure the interaction potential, Wint(?') 
with r = X2 — xi, between a pair of functionalized parti- 
cles sitting at positions (xi, X2) in a bias potential (laser trap) 
of strength Ubms{xi,X2)- We achieve this goal by performing 
two experiments (Fig. 2a). In one experiment we place two 
colloidal stars in two separate laser traps and in the other ex- 
periment we place two bare colloids in the same two traps. For 
both experiments we measure the separation histogram of the 
colloids as a function of the trap separation. The potential of 
mean force, M^sub, is then obtained by subtracting the results 
from each experiment. 

Wsn^if)/kBT = - log[/f(f)] + l0g[/nf(f)] (1) 

with /f(r) and fnf{f ) the fraction of measured displacements 
that fall within the histogram bin associated with the displace- 
ment value r for functionalized and non-functionalized beads, 
respectively. While this subtraction method has been used in 
previous experiments lfTll[T2ll . we rigorously prove its validity 
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FIG. 2: (color online). Excluded volume interaction of anchored 
rods, (a) the schematic and (b) the fluorescence image of phage- 
grafted beads in optical traps, (c) separation histograms of (A) bare 
beads and (B) phage-grafted beads for the same trap locations. The 
scale bar in (b) is 1 nm. 



here and show how to implement it over multiple windows. 

The fractions of measured displacements are governed by 
the Boltzmann distribution and given by 

/„f(r) = Z-,' J dx, J dx2e-^--(^^'^^)/'=»^^(xi - X2 - f) 

(2) 

and 

/f(f) = Zf^ J dxi J (1X26-^^'-'^''^''='^/''^'^ 

xe-i^.n.(=«2-:.i)/fcBTj^^^ - X2 - f) (3) 

with 6{r) the Dirac delta function and the normaUzation fac- 
tors are given by 

Z,rf = J dxi J da;2e-^'"-(^i'^=)/'=«^ 



We change the integration variables to xi and r = X2 — xi 
and integrate over r to obtain 



Inserting this result into Eq.[T]gives the calculated potential of 
mean force: 

VK,ub(r) = VK,„t(r) + fceT log(Zf/Z„f). (6) 

We see that Wsub{f) = Win({f) plus a constant. As dis- 
cussed above, the strength of the laser traps, C/bias(a;i, X2), is 
such that the colloids sample only a small range and therefore 
only a small piece of the interaction potential Wint(r) is ob- 
tained. To determine a wider range of Wint(r) the laser trap 



separation is varied and Wsub{f) is obtained anew. Although 
the constant term is different for each separation of the traps, 
the entire potential can be stitched together to within a single 
additive constant by assuming that Wint{f) is continuous. 

As a check of this implementation of the umbrella sampling 
algorithm, we used computer simulations to model the experi- 
ment. The results validating this method are shown in Fig.|3}l. 

The experimental system is shown schematically in Fig.|2^. 
The fluorescence image of trapped beads is shown in Fig.|2]3. 
Optical tweezer setup is built around the inverted fluorescence 
microscope. A single laser beam is time-shared between two 
points via a pair of orthogonally oriented paratellurite (Te02) 
acousto-optic deflectors (AOD, Intra-Action, Bellwood, IL). 
About 30 mW of a 1064-nm laser (Laser Quantum, Cheshire, 
UK) is projected onto the back focal plane of an oil-immersion 
objective (lOOx, N.A.= 1.3, Nikon) and subsequently focused 
into the sample chamber. Spheres are trapped 5 pim away from 
the surface to minimize possible wall effects. We choose a 
set of umbrella window potentials by systematically varying 
the locations of the traps' centers rd and rc2- For each win- 
dow potential, six minutes of video are recorded for a pair 
of phage-grafted beads, and the separation probability distri- 
bution, /f(r) is obtained. It is a simple Gaussian if the sep- 
aration is large and the beads are not interacting. The dis- 
tances between trap positions are selected so that there are 
sufficient overlaps between adjacent positions. We collected 
data for ~ 30 different trap positions with 50 nm increments 
in separation to cover a wide range of the interparticle po- 
tential. Under identical conditions (microscope illumination, 
laser power, sample buffer, etc.), the experiment was repeated 
immediately for a pair of streptavidin-coated PS beads with- 
out attached virus to measure fnf{f). For all experiments, sta- 
tistically independent configurations of beads were sampled at 
30 frames/sec. We analyzed the video images using a custom 
program written in the language IDL 1 16]. By constructing a 
histogram of center-center separations on 10"' images in each 
window, we found clear differences between the separation 
probability distributions of virus-grafted beads /f (r) and bare 
beads /„f(f) (Fig.El;). 

Fig. I3J) shows the interaction potentials measured between 
two M13-grafted microspheres with varying solution ionic 
strengths. The interactions are all purely repulsive. The po- 
tential decays to zero as the distance between sphere surfaces 
increases beyond twice the virus length. There is a strong de- 
pendence of the pair potential on the ionic strength of the sur- 
rounding medium. A decrease in the solution ionic strength 
leads to increased interaction between spheres grafted with 
charged rods. We compare the interaction potential between 
microspheres at grafting densities of 80 and 135 viruses per 
sphere (Fig. The increase in density by 68% increases 
the pair-potential by a factor of 2.6, but does not change its 
functional form. 

We calculated the interaction potential arising from the os- 
motic pressure of counterions trapped within the grafted lay- 
ers based on the mean field calculation theory of Jusufi |6, 7], 
except modified for the case where the density of fixed charges 
on the grafted rods is small compared to the salt concentration. 
In particular, the densities of positive and negative ions within 
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FIG. 3: (color online), (a) A set of interaction potentials of M13-grafted micron-sized polystyrene spheres acquired from each umbrella 
window, with the ionic strength 1—14 mM and the grafting density a — 135 phages/bead. The solid symbols indicate the potential extracted 
from the histograms shown in Fig. (b) Pair interaction potentials of colloidal stars at varying solution ionic strengths with a — 135 
phages/bead. Symbols: experiment; dashed lines: theory; solid lines: single exponential fits. (o,x') and (x): 2.8 mM; (A,y') and (y): 14 mM; 
(□,z') and (z): 28 mM. (c) Pair potentials at different grafting densities with 1—14 mM. (o,x') and (x): 135 phages/bead; (A,y') and (y): 
80 phages/bead. (d) Interaction potentials U/knT = Be^^"'' employed in the Langevin dynamics simulation (Solid lines) and bias 
potential L'bias(i'i, 1*2 ) = |fci|ri — rci|'^ + |fc2|r2 — rc2f- Pair potentials extracted using the umbrella sampling (empty symbols). B = 6 
(circle) and 20 (square). Insets: data replotted to facilitate comparison. ^ is the separation between the surfaces of spheres and L the virus 
length. 



the grafted layer p±(r) are given by p± « ps ± 0.5/9f(r), with 
Ps the salt concentration and pf{r) = XfNf/ {Airr'^) the fixed 
concentration of negative charges on the grafted rods, with 
N{ the number of rods per colloid and Af = 1.7e^/nm an 
adjustable parameter for the linear charge density renormal- 
ized by condensation. The counterion excess free energy is 
calculated by integrating over the volume of the grafted layer 

dOsmeJI^^^^''^^'^'^ drr'^Pf{rf/{2p,) with R, = 0.5 
the core radius, 9 the angle with the center to center vector 
for the pair of colloids, and = C + 2i?c the center to cen- 
ter distance. Following the interaction geometry depicted in 
Fig. 5 of Ref. |(7J, the height of the grafted layer is Lg = L 
for e > 00 and Lg = {2 cos 9) - for 9 < 9o with 
cos 6*0 = 0.5^c/(i + ^c)- We have assumed pf <C Ps and 
that the height of the unperturbed grafted layer is equal to the 
length of a virus, L = 880 nm (i.e. we neglect rod orien- 
tational fluctuations). Theory and experiment are compared 
in Figs. [3]5 & c; the agreement is rather good considering 
five measured potentials are fit with one value of the effective 
charge. 

We also calculated the interaction due to the excluded vol- 
ume of grafted rods based on the Onsager second virial expan- 
sion of the free energy [17J. For two rods, each with a spec- 
ified orientation, we find the pairwise excluded area, or the 



space of relative grafting locations for which the rods overlap. 
The calculated interaction potential due to rod excluded vol- 
ume was significantly smaller compared to the interaction due 
to counterion osmotic pressure; the latter interaction agrees 
well with the measured interaction potentials. 

In conclusion, the umbrella sampling method is applied 
to extract the pair potential of the colloidal stars which are 
trapped with optical tweezers. The method allows for mea- 
surement of potentials of the order of 100 k^T, an energy 
much greater than previously measured with line traps. The 
large measured repulsive energy between colloidal stars is 
consistent with the osmotic pressure of counter-ions between 
the charged brush, while a second virial theory based on the 
Onsager approximation significantly underestimates the pair 
potential. The construction of colloidal star polymers from 
genetically engineered viruses opens the possibility of a sys- 
tematic study of hybrid colloidal materials exhibiting complex 
phase behaviors. 
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